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Abstract 

The mass transport by a Burgers velocity field is investigated in the framework 
of the theory of stochastic processes. Much attention is devoted to the limit of 
vanishing viscosity (inviscid limit) describing the "adhesion model" for the early 
stage of the evolution of the Universe. In particular the mathematical foundations 
for the ansatz currently used in the literature to compute the mass distribution in 
the inviscid limit are provided. 
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1 Introduction 



It was proposed by Zeldovich [|TJ, that a possible model for the description of large scale 
dynamics of the mass distribution, in the early stage of evolution of the Universe, is 
provided by: 

3(V + v-Vv = fl, v (x,t = 0) = Vy(x) 

d t p + div(vp) = , p(x, t = 0) = p (x) U 

Here, the mass density is driven by a rotation free velocity field, solution of the Riemann 
equation. 

As it is well known, the Riemann equation with a gradient initial condition has multi- 
stream solutions H,@. This phenomenon is easily understood in the Lagrangian picture 
of the flow: the point-mass particles evolving along the characteristics, collide after finite 
time. The study of the multi-stream solutions requires a very subtle analysis of the 
caustics formed after a finite time by the Lagrangian solutions of the Riemann equation 

!• 

One way to avoid multistreaming is to introduce the "adhesion model" (see j3], ||, 
0, [0 with references therein). 

a (V + vVv = i/Av, v(x,i = 0) = W(x) () 
d t p + div(vp) = , p(x, t = 0) = po(x) U 

where the velocity field is now governed by the Burgers equation ||. 

The introduction of a small diffusion term is expected to have a smoothing effect only 
in the regions where shocks are about to occur: the limit of vanishing viscosity of the 
Burgers equation, the inviscid limit, selects one solution for the Riemann equation. 

The main disadvantage of the adhesion model (0) is the loss of a unique Lagrangian 
picture for both the velocity field and the mass density. The consequences are not only 
conceptual but also practical, since for (|2]) the solution of the Burgers equation does 
not provide immediately a definite algorithm which solves the continuity equation. The 
problem is therefore moved to the mass density evolution. 

In the one dimensional case, because of the strong topological constraint, there is 
general agreement on the idea that for times long enough that shocks appear, an initial 
uniform density field evolves into a singular distribution which describes the formation 
of point-masses on a background of a smooth, diluted, density field. The point-masses 
are situated at the shock positions of the solution to the Burgers equation and may be 
referred to as macro-particles. The mass of a macro-particle is equal to the integral of 
the initial mass density extended over the interval of initial positions which fall into the 
given shock at some time prior to that under consideration. 

In the two dimensional case the situation seems to be more subtle. Here we have the 
appearance of ruled surfaces, ribbons, where tangent planes touch the graph of the convex 
hull at a segment instead of a point and of triangles, corresponding to tangent planes 
with triple contact. In the Eulerian plane to these regions are respectively associated 
shock-lines and shock nodes. 
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A further source of complication is that a physically reasonable || initial condition 
for the Zeldovich model (Q) is to assume the velocity potential, at time t = 0, given by a 
generalized Brownian motion in the space variables. 

According to the adhesion model (§), in order to compute the mass density distribution 
at any t, first, we should integrate the rotation free Burgers equation and then use the 
resulting velocity v = Vfi field to solve the ordinary differential equation (ODE) 

x t = Vfi(x t ,t; v) , x = a (3) 

as a function of time and of the initial condition. Then the limit of vanishing viscos- 
ity should be taken to define the inviscid (weak) limit of every mass averaged smooth 
observable. 

It is clear that such a program in the realistic three dimensional case is very resource 
demanding for its numerical implementation. Therefore the question naturally arises 
whether a regularized Lagrangian picture exists which simultaneously solves both the 
velocity field and the mass density equation. 

In 0, Vergassola et al. proposed a numerical algorithm (the VDFN algorithm from 
now) based on the generalized Legendre transform of the potential field Q(pc, t;v = 0). 
The underlying ansatz is that the values of x, performing the transform, define a solution 
of®. 

In recent works |§, [[?]] it has been proposed that a consistent regularization of ([[]) 
with a simple numerical implementation, can be found by introducing a small viscosity 
term also in the right hand side of the continuity equation. 

In the present paper it is shown that a Lagrangian picture is uniquely associated to 
the Burgers equation. The characteristics of the Burgers equation are the realizations of 
the solution of a stochastic differential equation (SDE) describing a backward diffusion. 

In such a context, the mass transport along the characteristics of the Burgers equation 
is described by a backward Fokker-Planck equation, whose solution is related to a difficult 
inverse problem when the given boundary conditions are the initial velocity potential and 
the initial mass density distribution, (section ||) . 

On the other hand the Lagrangian approach shows that the choice of the Burgers 
equation in order to regularize the Riemann equation, is not the natural one and it is 
"exact" only in absence of shocks. In section |3] the natural regularization is introduced. 
It is shown that it provides a simple algorithm in the inviscid limit for the solution of 
both the velocity field and the mass density. The limitation of such a procedure is that 
it imposes, in order to be exact, some restrictions on the initial conditions of the velocity 
field. 

In section [| the approaches previously introduced in the literature (0,0) are reviewed 
in light of the results provided by stochastic calculus. In particular it is shown that the 
VDFN algorithm corresponds to an "effective" mass transport along the trajectories of 
the backward diffusion underlying the Burgers equation. 

Furthermore the pair of PDEs specified by the Burgers equation together with the "ef- 
fective" mass transport equation, has the Zeldovich model as the weak limit for vanishing 
viscosity for any initial data. The price to pay is in the artificial nature of the procedure. 
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The results of the present analysis suggest that both from the physical and the al- 
gorithmic point of view, the more natural way to introduce adhesion models is that one 
based on conservation laws [|T0 . 



2 Lagrangian picture for the Burgers equation 



Let us consider the backward stochastic differential equation ||11| , fl^, []13|l ) |T4 

= — v(u) s , t — s) ds + \[2v dw s s < t 



d s uj s 

U = X 

The Ito differential of the drift field along the flow is: 
d s v(u) s ,t - s) = 

= [— d t ~ s v(u} s , t — s) — v(o; s , t — s) ■ Vv(w s , t — s) + z/Av(a; s , t — s)]ds + 
+ V2udw s ■ Vv(u s , t — s) 



(4) 



(5) 



where dw t defines, as usual, the stochastic differential of Brownian motion. 
If v(x, t) satisfies the Burgers equation, the equality holds: 

v(x,t) =< v(^ x '°,0) > 



(6) 



where average < . . . > is taken over the realizations of the diffusion defined by (|J). 

The physical meaning is clear: the Burgers equation defines a velocity field constant 
on the average over the random trajectories of (f|). The conservation law allows a straight- 
forward integration of the rotation free Burgers equation. 

The basic object to be considered [15| is the transition probability defined by: 

P(_)(y,s | x,t) =< S(u - x) 8(u t ^ s - y) > s<t (7) 

The transition probability (|?|) satisfies [15| in the variables (y, s) the backward Fokker- 
Planck equation: 

d s p { -) + div(vp ( _)) + ^Ap(_) = (8) 

with the final condition: 

limp ( _)(y, s|x,t) =<y(y-x) (9) 
Let us introduce the auxiliary stochastic process described by the SDE: 



dz s 

Z() 



(10) 



The transition probability density for this process at time s is easily found to be Gaussian 
with expectation value x and variance 2 v s. 
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We can transform (0) into an average over the realizations of z f : 

P(-)(y> s I x ^) =< <K z o - x) 5(zt_ s - y) (zts) > (11) 

The Jacobian of this transformation is supplied by the famous Girsanov formula (see for 
example ]T3| or fT6|j ): 

[xt-s) = exp[- / t= dw u - / rfwj (12) 



dP Zf _« Jo \/2v Jo 



where the relation v = Vfi has been used. Furthermore along the trajectories z s the Ito 
differential of the velocity potential becomes: 

d s Q(z s , t — s) = [— d t - s £l(z s , t — s) + v AQ(z s , t — s)] ds + V2v dw s ■ Vf2(z s , t — s) (13) 

If we use to eliminate the stochastic integral in flT2"D and we impose the normalization 
condition ^ 

/ d D yp^ ) (y,s\ X ,t) = l (14) 

J — oo 

we find the velocity potential to be 

fi(x,t) = -2Wn{(--i-)T r d n yeM Jj^Z^l_YM ]} (is) 

and 

exp r (y-x) 2 _ fXMl 

/ I + \ P[ 4i/(t-s) 2!/ J / 1c x 

p ( _)(y,a|x,t) = — -— )2 ^y- (16) 



It is easy to verify that ( |T6"D is a well defined Markovian transition probability for any 
s < t. 

In the physical literature is more common to deal with stochastic calculus in the 
Stratonovich representation which has the advantage to preserve the ordinary rules of 
differential calculus ]17| . 

If we express the stochastic Ito integral, appearing in (0), in terms of the correspond- 
ing Stratonovich's [Tj|, then what we have done is nothing else than proving that the 
transition probability (|7j) is given by a Feynman path integral with Lagrangian: 

r/ • \ ||X+ Vfi(x,t - S)|| 2 1 A _,. 

4^ 2 



From (|17f) the derivation of (|Tq) is then trivial. 

The solution of the Burgers equation at arbitrary time takes the form: 



/oo 
d D j/W(y)p(_)(y,0|x,t) 
-oo 
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The knowledge of (|7|) solves also the problem of the forward evolution of a passive 
scalar driven by the Burgers equation. More explicitly if we consider: 



a t C + v-VC = vAC + E 

C(x,0) = C (x) (19) 

where E = E(pc,t) is an external forcing, then we have 

/oo ft roo 

d D yCo(y)p i -)(y,0\x,t)+ / / d D y du E(y, u)p H (y, u | x, t) (20) 
-oo JO J —oo 

As a function of (y, s) the transition probability (|7|) describes the mass transport along 
the characteristics (|4]) of the Burgers equation. 
The equation with final condition 

d s p + div(vp) + vAp = , p(x, t = T)= pr(x) (21) 

is readily solved for any t < T by: 

/oo 
d D zp ( _)(x,t|z,r)pr(z) (22) 
-oo 

On the other hand the boundary conditions for the Zeldovich model provide the initial 
mass distribution. In order to solve the mass density evolution in the interval [0, T], we 
have, first, to solve the inverse problem: 

, N r Wx). r°° d D z . (x-z) 2 fi(z,T). . 

Po (x)exp[^]= / exp[-i-— L + _L^i pZ)T (23) 

In principle, once we have solved equation (p3|), it would be possible to specify the 
mass density evolution in the interval [0, T\. The fact that the interval is finite, is not, 
from this point of view, a limitation since p(x, T) could be consistently used as a new 
initial condition at time t — T to iterate the procedure, as far as the integral in ( p3|) is 
convergent. 

It is worth to note |13 |, [14|, that the problem (p3|) can be reformulated by looking 
for the conditions that insure the existence of the forward stochastic process specified for 
any t G [0,T] by the probability density, solution of (pi]) with p(x, T), the unknown final 
condition to be determined. 

The crucial point is that, in general, the transition probability density (0) does not 
have an inverse: because of the non-reversible nature of the diffusion dynamics, it satisfies 
only a semi-group property [ 15 ]. 

On the other hand, it is interesting to understand the physical meaning for the occur- 
rence of the inverse problem (P3"D, when our starting point is to find a coherent Lagrangian 
regularization for the Zeldovich model. 
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3 The forward diffusion 

As it is clear from the above discussion, Burgers equation is equivalent to the SDE: 

d s u s = -v(u 8 ,t- s)ds + V2udw s , u = x , . 

d s v(u s , t — s) = ^/2vdw s • Vv(u s , t — s) , v(x, 0) = v (x) 



(25) 



The limit of vanishing viscosity of (l24f) is 

<i s u; s = — v(u; s , t — s) ds , uj = x 

d s v(uj s ,t - s) = , v(x,0)=v (x) 

On the other hand, the characteristic equation of the Zeldovich model (Jl]) is 

d t £t = v(£t,t)dt, £o = y / 9fi x 

dtv(6,*)=0 v(a,0)=v (y) W 

If £ t y '° is a solution of fl26|) in [0, i] then dJ s x, ° = £ t ^* is a well defined solution of d25|) on/?/ 
provided that the solution for (|26|) exists and is unique in [0,t\. 

If the latter condition is satisfied, (|25|) , allows to solve the inverse problem implicit in 
the Riemann equation or in the equations for any passive scalar conserved along the flow: 

dMHt y ' ,t) = 0^v(tf> ,t)=v(y,0) (27) 

Namely, if v is well defined for all s G [0, t] the Riemann equation is equivalent to 

d a v(x,t- s) - v(x,t- s) • Vv(x,t- s) = 0, v(x,0)=v (x) (28) 

The straightforward consequence is the well known fact [@ that for any t such that no 
collision occurs between the Lagrangian particles, we have: 

v(x,t) = Vo(^ y '°' ) 

= x-v (y)t (29) 

It is worth to note that and the steps leading to it are exactly the deterministic 
counterpart of (§) and of the procedure providing the integration of the Burgers equation 
in the particular case of rotation free initial conditions. 

Beside the velocity field, the Zeldovich model requires the mass density field. It is clear 
already from the deterministic case that the knowledge of u; t x '°, alone, does not provide 
the solution of the continuity equation: 

d t p + div(vp) = , p(x, t = 0) = p (x) (30) 

The only expression consistent with ( |3~0"D we can construct, is 

/oo 
d D x6(y-u t y 2° s )p(x,t) (31) 
-co 
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where p(x, t) is a given solution of the continuity equation for t > s. Again the determin- 
istic case is a mirror of the situation we faced studying the mass transport by the Burgers 
equation, i.e. equation (f?3|). 

The choice to regularize the Riemann equation by resorting to the Burgers equation 
is equivalent to select a solution for the Riemann's, when shocks occur, by extending 
the time reversal conjugation between the exact equation for the characteristics (^) and 
the deterministic limit of (p4"|). This procedure is of clear advantage when dealing with 
the velocity field alone, but the price to pay consists in the difficulties which arise when 
turning to the problem of the mass evolution. 

The alternative approach is to try to regularize the Riemann equation by looking to 
the stochastic generalization of its Lagrangian picture. 

The SDE to be considered is then: 



d& = v(£ tj t) dt + V2udw t , £,o = y /o 2 \ 

ci t v(6,t) = v / 2^(iw t -Vvfe,t), v(x,0)=v (x) 1 ] 



Hence, the system of PDEs implied by (|32|) is 



d t v + v- Vv + i/Av = 0, v(x,t = 0) = W(x) 

d t p + div(vp) = z/Ap , p(x, t = 0) = p (x) 

For a rotation free initial velocity field, the integration of the pair of equations ( p3|) is 
achieved in terms of the transition probability density defined by ([32|), following the same 
steps leading to (|16|). 

The result is, for any s < t: 

exp[-i^i + ^l 

/ I \ ^ L iu(t-s) 2 v J / q q \ 

P(+) (x, 1 1 y , s) = — — 2 5( - (33) 

where the velocity potential is fixed by the normalization condition and is equal to 

1 ,d f°° n (x-z) 2 S(Z,T),- . , 

S(x,t =2flog{b 77^ tt] ~ / rf D zexp- . 1 , T ; ,x +-V-^> 34 

4tiu(1 — t) J-oo 4is(l —t) Iv 

Therefore, in order to achieve the solution in the interval [0, T] we need to solve the inverse 
problem: 

T V(x), f 00 d D z r (x-z) 2 S(z,T), , s 

exp — ^ = / 5-exp - - J + V ' ; 35 

PL 2// J 7-00 (47rz/T)f PL 4z/T 2u 1 K ' 

Equation (|35| ) is the natural stochastic generalization of the inverse problem fl2"T|). Once 
5(x, T) is known, it can also be assumed as the new initial condition at time T to iterate 
the procedure. 

In the limit of vanishing viscosity, equation (|35|) takes the form: 

V(x) = sup[H(z,T)-^^] (36) 
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The latter expression is basically a (generalized) Legendre transform |J: 

T V(x) - ^ = sup[TH(z, T) - y + x • z] (37) 
It is tempting to solve (|36|) as: 

H'(x,T) = inf[y(z) + ^£] (38) 

This solution is not exact for a general initial condition: due to the convexity properties 
of the Legendre transform, when we substitute S(x, T) back in (|36| ) what we find is not 
V^(x) but its convex hull. Only for this latter the transform (|3l)D is involutive. Therefore, 
by assuming S(x, T) = H'(x, T) we implicitly adopt the consistent initial condition. 

If we accept, despite of the above restrictions, the pair (^) as a regularization for 
the Zeldovich model then we are provided a simple algorithm for the computation in the 
inviscid limit of both the velocity field, by means of fl34]) , and the mass distribution, by 
means of: _ 

/oo 
d D y P(+) (x,t|y,0)p (y) (39) 
-oo 



4 Non Lagrangian approaches 

In M it has been proposed to regularize the Zeldovich model with the pair of PDEs: 

a t v + v • Vv = z/Av , v(x, t = 0) = W(x) 

+ div(vp) = /iAp , p(x, t = 0) = p (x) 

The meaning of ( |40D is to assume the velocity field solution of the Burgers equation as 
an external field in the Langevin equation describing the motion of a point mass particle: 

d t y. t = v(ya t ,t) dt + y^dwt , x = y (41) 

Hence, the acceleration felt by the point-mass particle is 

dtv(pct,t) = (Li + v)Av(x t ,t)dt + ^T/ldwf Vv(x t ,£) (42) 

It is worth to note that fl42|) predicts for the case /j, = v a a mass dynamics different 
from ( |32"D even in the limit of vanishing viscosity: only in absence of shocks the two 
dynamics become equivalent. 

In the general case the solution of mass distribution resulting from ( fTO|) is difficult. 
The path integral approach clearly shows that it is equivalent to an Euclidean Schrodinger 
equation in a potential given by the Laplacian of the velocity potential of the Burgers 
equation. Namely, the transition probability for the mass density in (|40|) is 

Q(x,t) n( y ,o) 

p(x,t|y,0) = e^K(x,t; y,0)e~5T (43) 
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where K is given by 



#(x,t;y,0) = P X Vx u exp{- f + ^-^AQ(x u ,u)) du] (44) 

ix =y JO 4 fj, 2 jj, 

and it satisfies 



<9 t K + X Afi = /iAK (45) 

2 // 

hm^(x,t;y,0) = <5(x-y) (46) 

There are two evident cases when the path integral (|44]) is of practical use. When the 
Laplacian of the velocity potential does not depend on the trajectory, all the integrations 
turn out to be Gaussian. This condition is satisfied for: 

1. V(x) = vo • x. Then the transition probability is: 

(x-y-vp t) 2 -\ 



cxp [ 



which, for every smooth enough initial mass distribution, gives for fj, J, 0: 

p(x,t) = p (x-v t) (48) 
in accordance to Galilean covariance. 
2. V(x) = J^eL. By means of (33) we get into 

In the inviscid limit for smooth enough distribution we obtain: 

The peculiarity of these two examples appears also in the fact that the results do not 
depend on the viscosity v of the driving velocity field. Actually the initial conditions for 
the velocity potential given above, are such that no shock appears at any time and the 
Laplacian of the velocity field is zero. 

More results can be derived for the case fj, — v. The crucial observation is that the 
fundamental solution of the mass density problem can be rewritten as 

dtp + v • Vp = vA p — pV ■ v 
limp(x,t|y,s) = £(x-y) (51) 

tls 

This means that the transition probability p can be computed as an average over the 
random trajectories, solutions of (|4j) 

p(x,*|y,a) =< J(M, M) <K"i-° - y) > s<t (52) 
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where 

rt-s 



Jit, s, {u}) = expf- / S (JmV-v«'V-«)] (53) 
Jo 

is the Jacobian of the change of variables in the functional integration between the solu- 
tions of the backward diffusion fl24[) and those of ( f4~l~D . 

In one dimension this observation leads to a straightforward integration of (|5T]) : 



< J(t, s, {u})5(u t x J° s - y) >=< - y) d x u^ s >= -d x < 6(y - u*§ > (54) 

where 9 is the step function. It follows: 

ry 

p(x,t\y,s) = -d x dzp(-)(z,s\x,t) (55) 

J — oo 



One can easily check that (55) verifies ([H]) and it is a Markovian transition probability. 

Unfortunately, the identity in ([54]) does not hold in more than one dimension. 

For the general case, in J7| it has been introduced a "mean field approximation" for 
(pip. Here a different interpretation of the result is proposed. 

The mean field theory of M is equivalent to substitute to (pTj) the equation 

d t p* + v • Vjt>* = v A p* — p* < V • v > 
limp*(x,t|y,s) = 5(x-y) (56) 

tls 

where the average means 

/oo 
d D xV ■ v(x,t)p*(x,t|y,s) (57) 
-oo 

Furthermore we need to impose, for every t < s, the constraint: 

/oo 
d D xp*(x,t \y,t) = 1 (58) 
-oo 



The integration of ( |56| ) is immediate. Since (|5T| ) does not depends on x, the corre- 
sponding term can be extracted from the path integral, which is reduced to an average 
over the realizations of (P4"|): 

(x, 1 1 y , s) = e~ /. du < v ' v> p(_) (y , s | x, £) (59) 
The normalization condition then fixes the value of the prefactor. Finally we get into: 

expf (x ' y)2 + 

p*(x, t y,s) = — ^ — Ta n , n (60) 

J-oo" z e3C PL 4i/(t-«) ^ 2i/ J 

Here the velocity potential of the Burgers equation fi(x, i) explicitly appears. The 
effect of the average (|57]) is to define an "effective" mass transport along the characteristics 
of the backward diffusion (124]). The "effective" theory becomes an exact solution of the 
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continuity equation in the (weak) inviscid limit: for any smooth observable described by 
a scalar function /(x), it is readily verified that 

/oo 
d D x d D y /(x)p„(x, t\y,s\ v)po(y) (61) 
-oo 

satisfies in the weak sense the continuity equation and we have, 

/oo 
d D yf[x(y,t)}p (y) 
-oo 

x(y,t) = ar gsup[fi(z,t;z/ = 0)-^f^] (62) 
z Z T 

where fl(x, £ ; v — 0) is the inviscid limit of velocity potential of the Burgers equation. 

The latter result ( p2|) is exactly the ansatz of the VDFN algorithm used in , in order 
to compute the evolution of the mass density field driven by the Burgers equation. It is 
also worth to note, that in the case of a convex differentiable velocity potential, equation 
(B3) implies free motion for the point-mass particles. 



5 Conclusion 

The "natural" regularization ( |33|) of the Zeldovich model by means of the introduction 
of a small viscosity coefficient leads to the inverse problem fl35|) which is the direct gen- 
eralization of the one occurring in the deterministic case. The inviscid limit drastically 
simplifies the situation, although it imposes some restrictions on the initial conditions. 
Nevertheless, if we neglect such difficulty, the forward diffusion approach (|33D provides 
us a simple algorithm to compute, in the inviscid limit, the velocity field and the mass 



distribution at any time (equations (j34]) , (j38D and (|39|)). This procedure is exact for convex 
initial conditions. 

To select the solution of Riemann equation by extending at larger times the corre- 
spondence with the backward characteristics of the Burgers is of real advantage only if 
we are interested in the velocity field alone. 

The occurrence of inverse problems can be avoided if we associate to the Burgers 
equation the "effective" equation (EH) for the mass transport along the trajectories of the 



backward diffusion (p4|) . The solutions (|T5| ) and (^) of the pair of equations specified 



by the Burgers together with its "effective" mass transport tend in the limit of vanishing 
viscosity to a weak solution of the Zeldovich model. 

For the mass density field, such solution (equation (^)) was proposed as an ansatz in 
H] . The algorithm defined by (^) has the advantage of an easy numerical implementation 
on a computer like that one provided by the "natural" regularization (|33|) . Furthermore 
in comparison with the latter it does not impose any restriction on the initial conditions. 

The main disadvantage of (|62|) is its intrinsic non-locality which makes its use artificial 
from the microscopic point of view. 

The interpretation of the transition probability (|60|) , defined by the "effective" theory 
(pBj) , as an approximate solution of the exact mass transport (^Tj) by a forward diffusion 
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with external drift field given by the Burgers equation, is certainly correct only before the 
occurrence of shocks. For larger times the characteristic equation (fy]) with /x = v seems 
to indicate a different behavior over shock domains. Only in one dimension, heuristic 
arguments can be provided, 0, to show the equivalence, in the weak inviscid limit, of 
(|BT|) with the mass transport described by the exact transition probability ([55]). 

Finally it must be remarked that starting from the basic kinetic equations, the Zel- 
dovich equations (|1|) are not the unique starting point for the construction of adhesion 
models. Let us consider 

dtf + p- V x / = (63) 

where / = /(x, p, t). The ansatz 



/(x,p,t)=p(x,t)%-v(x,t)] (64) 



leads to the equations 



pd t v + p (v • V)v = 
d tP + V ■ (pv) = (65) 

By means of the continuity equation, the system can be recasted in the form: 

d t {pvi) + Vj(vj pvi) = 

d t p + V J (pv j ) = (66) 



In [[L(| , it is proven that the latter pair of PDEs is equivalent to the Zeldovich model 
only in absence of shocks or for an uniform initial mass distribution. 

The simple meaning in terms of conservation laws of equations (^B|) allows to achieve an 
exact algorithm for the solution which turns out to be of simple numerical implementation 

0- 
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